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Within the framework of perturbative resummation scheme of Pis- 
^ ■ arski and Braaten, the decay- or damping-rate of a moving heavy quark 



(muon) to leading order in weak coupling in hot QCD (QED) is exam- 
ined. Although, as is well known, the conventionally-defined damping rate 
diverges logarithmically at the infrared limit, shown is that no such diver- 
gence appears in the physically measurable decay rate. The cancellation 
occurs between the contribution from the "real" decay diagram and the 
contribution from the diagrams with "thermal radiative correction" . 
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In hot gauge theories, a consistent perturbation scheme (the hard-thermal-loop 
(HTL) resummation scheme) has been established [fll £3]. This scheme enables con- 
sistent evaluations of any thermal reaction rates to leading order in the coupling 
constant. The HTL resummed propagators screen the infrared singularities, and ren- 
der otherwise divergent physical quantities finite, if they are not sensitive to a further 
resummation of the corrections of 0(g 2 T). Much interest has been taken [p], |3], || in 
the so-called damping rate of a particle moving in a hot quark-gluon (electron-photon) 
plasma, because this is an example of those which are sensitive to the further resum- 
mation mentioned above. In other words, the damping rate of a moving particle is 
still divergent even within the HTL resummation scheme. It is anticipated |1|, f|] that 
the further resummation makes the damping rate finite. 

In this Letter, within the HTL resummation scheme |I|, |2J, we analyze the decay- 
or damping-rate of a moving heavy quark (muon) in a hot quark-gluon (electron- 
photon) plasma, in terms of a measurable quantity. We shall find that, although the 
damping rate as defined in a conventional way diverges, no such divergence occurs in 
the physically measurable decay rate. The mechanism of resulting the finite decay 
rate is the same as the familiar Bloch-Nordsieck mechanism, which operates to cancel 
the infrared divergences in reaction rates in vacuum theories ||. 

We consider the heavy quark with velocity v = p/E (0(g) << v) injected into 
the quark-gluon plasma. By the term "heavy quark" we mean so heavy a quark 
that it is not thermalized, i.e., e~ E l T << 1. The heavy muon immersed in the hot 
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electron-photon plasma may be treated in a parallel manner (see below). We employ 
the real-time formalism of thermal field theory |J, which is formulated on the time 
path — oo — > +00 — > —00 — > —00 — 2T™ 1 , in a complex time plane. 

The decay rate of the heavy quark with any one of the color states is obtained as 

M 1 

n d = -i--Y.u s {P)^i{P)u s {p) ) (i) 

where £ 21 (P) with P = (P, p) is the (21)-component of the self-energy matrix 
Ejj(P) = 1,2) of the heavy quark. Here i and j designate the type of heavy 
quark-gluon vertex || §], the type-1 (type-2) vertex stands for the vertex of physi- 
cal or type-1 (thermal-ghost or type-2) fields. [The "Feynman rules" in the above- 
mentioned real-time formalism is equivalent to the circled diagram rules of Kobes and 
Semenoff J7|, provided that the physical (thermal-ghost) field is identified with the 
field of "uncircled" ("circled") type in J7j.] To one-loop order in the HTL resummation 
scheme, £ 2 i is diagrammed in Fig. 1 with i — 1 and j = 2; 

S 21 (P) = -ifCf J ^*A^(Q)7^ 21 (P')7,, (2) 

where C f = 4/3 and P' = P - Q. In (2), 

5 21 (P') = -2m 6(p' ) (7 • P' + M) 5{P' 2 - M 2 ) , (3) 

and *A%i(Q) is the (21)-component of the (HTL resummed) effective thermal gluon 
propagator. For the heavy muon in the hot electron-photon plasma, in (2) and in the 
following, g should read e and Cf = 1. 
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Throughout this Letter, we keep only the terms that yield the leading contribution 
of 0(g 2 In g), which, in fact, diverges logarithmically (cf. (8) and (9) below). It is 
well known []T], |^, || that the relevant terms come from the chromomagnetic part of 
*A£i(Q) in the soft Q = (q ,q) region, < O(gT); 



[*A%(Q)r S = - ~^q P q a ) [0(qo) + n B (\q 

1 



x 



Q 2 - n t (g + iqoe,q) 



c.c. 



(4) 



where 5^ u = (0, 1, 1, 1), q = q/g, q — |q|, and "c.c." stands for complex conjugate. 
In (4), n.B ( | go I) is the Bose distribution function and |9|, || 



n t (9o,9) = 2 m 



Qo\ , Qo(q -%) i Qo + q 



,q 



2q 3 



qo-q 



rrin 



l{3 + ^)(gT) 2 forQCD 
| (eT) 2 for QED , 

with Nf the number of quark flavors that constitutes the quark-gluon plasma. Here 
we introduce the spectral density 



1 



1 



7T Q 2 - Tl t {q + te, q) 

By noticing that the factor 9(qo) in (4) can be ignored, since ned^o!) — T/\qo\ > 
0(l/g) » 1, we can reduce (4) to 



rA£XQ)) mag ~ -27TZ - P* ~8™q p q a ) - p t (q , q) . 

qo 



(5) 



Substituting (5) for *A£(Q) in (2) and using (3) with 6(P' 2 - M 2 ) ~ 5(v • q 



go)/(2P), we obtain 



E 21 (P) ~ ^ 2 C f | / ^ 1 *(v • q - g ) {^7o - (P • q) (7 ■ q) " M} p t (g , g) . 

It is also well known [0, [| £|] that the dominant contribution to E 2 i(P) or TZd 
comes from the region v ■ q = q << vq, where 



p t {qo,q) 



3 m% 



ggo 



4 g 6 + (? m?.) g 2 

Using this approximation for p 4 (go, g), we obtain 



E 21 (P) ~ ^g 2 C^m^(i?7o-M) 



d 3 g 



g 



(6) 



(27r) 3 g6+ ^ m ^ (v-q) 2 ' 
This expression is valid for g < O(gT) and v ■ q << 1. The integral in (6) diverges 
|], |3], [| logarithmically at the infrared limit q — > 0; 



« 2 - r 



E 21 (P) ~ g*C f -(E>yo-M) 



arctan 



7T^ p 

fO(m T ) dq 
x / 

■ g 

f - (7 2 C) - (E l0 - M) 
An p 



vm 



T 



(7) 



Inserting (8) into (1), we obtain 



Tl d ~—g 2 C f vT 



(9) 



The damping rate, 7, is defined [flj |3|, |J to be i times the imaginary part of the 
pole with Repo > of the quasiparticle propagator {P ■ 7 — M — E^P)} -1 , where 



Ejr(P) = -S 22 (po(l + ie),p) -S 2 i(po(l + ie),p) (cf. Chapt. 3 of §) . From (8) and 
(12)-(14) below, we obtain, for p ~ i£, 

Im£ 21 (P) ~ -2/mS 22 (P) 

~ ^{E l0 -M)n d . 



Thus, within the HTL resummation scheme, the damping rate 7 of the moving heavy 
quark is logarithmically divergent, the well-known result which became to be an 
opening of a variety of work jlL |3], f|] under the name of "damping-rate puzzle" . 

In an experiment, we "detect" the quark with momentum p' = p — q (Fig. 1). 
Then, the independent variables are p' and Z = pp'. Instead of using these variables, 
we use q = |p — p'| and z = p ■ q; p'Z = p — qz and p' 2 = p 2 + q 2 — Ipqz ~ p 2 — Ipqz. 
As mentioned above, the dominant contribution to 1Zd comes from the small region 
z = qo/{vq) << 1. Then, we consider the differential decay rate dTZa/dq, where the 
integration over z has been carried out. As is seen in (7), at q ~ 0, TZd exhibits a 
dq/q spectrum, causing the logarithmic divergence. It should be recalled here that a 
detector has a finite resolution with a typical resolution Aq. Then, what we measure 



Then we see that, to the present approximation, 



7 = 5*-. 



(10) 



as dlZd/dq 



q<Aq 



in an experiment is 
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(11) 

We will see below that the diverging decay rate (11) is compensated by the analogous 
term in the "elastic process" . 

Now we consider the contribution of the diagrams with "HTL-resummed virtual" 
gluon, Fig. 1 with i — j — 1 and with i = j = 2; 

K« virtuar = -i^£fZ s (P)[E n (P) + E 22 (P)]u s (P) 

= 9 2 ^C f j ^j:Eus(P)il»S u (P') lu *A%(Q)}u s (P), (12) 



=1 s 



where 



s ^ p> ) = (-y- \^'^ +M f v = L 2 ) • 



p/2 _ M 2 _ H 

As in the case of TZd above, the magnetic parts of *A^(Q) (£ — 1, 2) in the region 
| go I << <? < O(gT) dominates the integral in (12), where 

(*A^T(Q)) mag = -{(*A^(g)) mag }* 

= (5^ -~5^~5 va q p q a ) 
l + n B (\qo\) 



x 



n B (\qo\) 



Q 2 - n t (g + iqoe, q) 



Q 2 - n t (g - iqoe,q) 



(13) 



Substituting (13) with (5) for *A%(Q) {I = 1, 2) in (12) and using S U (P') + S 22 (P') = 
-2ni (7 • P' + M) 5(P' 2 - M 2 ), we obtain 



S n (P) + S 22 (P)~-S 21 (P), 



(14) 



and then 



virtual" 



(15) 



It is to be noted that the physically measurable quantity "at q < Aq" . 
[(dTZ/dq) (Aq)] q<Aq , is § 



dlZ , . ' 



q<Aq 



+ 



1 'R "virtual" 

dK d 



dq 



(A?) 



(16) 



q<Aq 



The factor 1 on the r.h.s. of (16) is the zeroth-order contribution. Substituting (11) 
and (15) with (9) into (16), we see that the cancellation of divergences occurs between 
"virtual"- and "real" -contributions, and find 



dq 



(Ag) 



1 -—g 2 C f vT 



lq<Aq 

Thus the measurable quantity, although sensitive to the resolution Aq, is free 
from divergence. The mechanism of cancelling divergences in (16) is exactly the same 
as in the cancellation of infrared divergences in vacuum theory (Bloch-Nordsieck 
mechanism) ||. In some reaction rate in QED at T = 0, both the diagram with 
electron (soft) bremsstrahlung and the diagram with radiative correction diverge, but 
they cancel each other: The former is the counterpart to [{dlZd/dq) (Aq) } q <A q (Eq. 
(11)), while the latter is the counterpart to TZ^irtuai" (Eq. (15)). It is worth pointing 



out, in passing, that we obtain 



dq 



(Aq) 



+ I -7 A dq=l 

lq<Aq JA 1 d( l 
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as it should be. 

Now we are in a position to discuss the physical content of (16) in terms of the 
processes taking place in the quark-gluon plasma. For any given diagram in the 
real-time formalism, general rules of identifying the physical processes are available 
||: Given a (thermal) diagram, like Fig. 1, representing some thermal reaction 
rate, one may divide it into two parts; the one is the reaction's S'-matrix element 
in vacuum theory and the other is the S*-matrix element. The S- and S*-matrix 
elements represent the reactions between the considered particle(s) (the heavy quark 
in the present case) and the particles in the quark-gluon plasma. The type-1 (type-2) 
vertices in the thermal diagram go to the vertices in the S- (S*-) matrix element. A 
thermal propagator with momentum K from a type-1 vertex to a type-2 vertex is 
"cut". When k > (k < 0) the "cut" is the "final-state cut" ("initial-state cut"). 
For a thermal propagator from a type-2 vertex to a type-1 vertex, the opposite rules 
apply. For more details, we refer to ||. 

We note that *A^i in (2) and *AgJT in (12) may be expanded in powers of n t 
(cf. (4) and (13)), and thus TZd and H" V irtuai" turn out to be represented by infinite 
series; lZ d = X^o^d an d ^"virtual" — EH^t'irtMr • From each contribution thus 
obtained, we can identify the physical processes according to the rules outlined above. 
An example of the physical process that is involved in TZ d {TV~^i rtua r ) ls depicted in 
Fig. 2 with the final-state cut line Cd (C„); the left side part of the cut line represents 
the S'-matrix element in vacuum theory, while the right part represents the S^-matrix 
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element. The group of particles on top of Fig. 2 stands for spectators which are 
constituent particles (quarks, antiquaries and gluons) of the quark-gluon plasma. In 
Fig. 2, Q is soft ~ O(gT), while K u K 2 , Pi and P 2 are hard ~ 0(T). It is to be 
noted that, in Fig. 2 with the final-state cut line C v , the heavy quark in the S^-matrix 
element is simply a spectator particle. 

For Q ~ 0, P ~ P', and then the "heavy quark detector" cannot discriminate the 
final-state heavy quark (P f ) in Fig. 2 with the final-state cut line Cd and the one (P) 
in Fig. 2 with C v . Then, what we measure corresponds to the sum of Fig. 2 with 
C d and Fig. 2 with C V) as well as of all other diagrams included in 1Z d , Eq. (1), and 



We like to emphasize again that two heavy quarks, one with P' — P (q — 0) and 
another with P' (|p' — p| < Ag), cannot be discriminated experimentally. In other 
words, one cannot recognize the heavy quark with |p' — p| < Aq as a (thermal) decay 
product, which is responsible for diverging damping rate 7, Eq. (10). On the basis 
of this observation, we propose to introduce a "observable damping rate" defined as 



This 7 bs(Ag) is the damping rate that originates from the heavy-quark decay into 
physically distinguishable states. 

It is to be noted that, for a heavy quark at rest, the damping rate 7 is infrared-safe 



K« virtuar , Eq. (12), (cf. (16)). 



(cf. (10)) 
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and finite quantity of 0(g 2 T) ||T], [UJ, and then the difference between 7 and •y \ ys (Aq) 
is negligibly small. The reason for the infrared safety of 7 is traced back to the fact 
that the heavy quark at rest couples only to the chromoelectric parts of *A2i(Q) 
and *A^(Q) (£ = 1,2), which, in contrast to the chromomagnetic part, develops a 
(thermal) Debye mass and there emerges the energy gap between the single heavy- 
quark state and the states with a heavy quark plus a "HTL-resummed gluon" with 
Q ~0. 
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Figure captions 



Fig.l. A thermal self-energy diagram of a heavy quark in real-time thermal field 

theory. % and j (i, j — 1, 2) designate the type of vertex. The blob indicates 
the (HTL resummed) effective gluon propagator. 

Fig. 2. Typical processes taking place in the quark-gluon plasma. The diagram with 
the final-state cut line Cd (C v ) is the process that is included in TZa (^"virtual") 
or in Fig. 1 with % — 1 and j = 2 (i = 1 and j = 1). The left side of the 
final-state cut line represents the ^-matrix element, while the right side rep- 
resents the S^-matrix element in vacuum theory. 
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